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Finite version of Random Domino Automaton (FRDA) - recently proposed in [l[ as a toy model 
of earthquakes - is investigated. Respective set of equations describing stationary state of the 
FRDA is derived and compared with infinite case. It is shown that for the system of big size, these 
equations are coincident with RDA equations. We demonstrate a non-existence of exact equations 
for size N > 5 and propose appropriate approximations, the quality of which is studied in examples 
obtained within Markov chains framework. 

We derive several exact formulas describing properties of the automaton, including time aspects. 
In particular, a way to achieve a quasi-periodic like behaviour of RDA is presented. Thus, based on 
the same microscopic rule - which produces exponential and inverse-power like distributions - we 
extend applicability of the model to quasi-periodic phenomena. 

PACS numbers: 45.70.Ht (Avalanches), 02.50.Ga (Markov processes), 91.30.Px (Earthquakes) 

Keywords: stochastic cellular automata, avalanches, exact solutions, toy models of earthquakes, forest-fire 

models, Markov chains 



I. INTRODUCTION 

The Random Domino Automaton, proposed in is 
a stochastic cellular automaton with avalanches. It was 
introduced as a toy model of earthquakes, but can be 
also regarded as an substantial extension of 1-D forest- 
fire model proposed by Drossel and Schwabl 

The remarkable feature of the RDA is the explicit one- 
to-one relation between details of the dynamical rules 
of the automaton (represented by rebound parameters 
[lijv defined in cited article and also below) and the 
produced stationary distribution rii of clusters of size i, 
which implies distribution of avalanches. It is already 
shown how to reconstruct details of the "microscopic" 
dynamics from the observed "macroscopic" behaviour of 
the system 0, Q . 

As a field of application of RDA we studied a possi- 
bility of constructing the Ito equation from a given time 
series and - in a broader sense - applicability of Ito equa- 
tion as a model of natural phenomena. For RDA - which 
plays a role of a fully controlled stochastic natural phe- 
nomenon - the relevant Ito equation can be constructed 
in two ways: derived directly from equations and by his- 
togram method from generated time series. Then these 
two results are compared and investigated in [(| 0] ■ 

Note that the set of equations of the RDA in a spe- 
cial limit case reduces to the recurrence, which leads to 
known integer sequence - the Motzkin numbers, which 
establishes a new, remarkable link between the combina- 
torial object and the stochastic cellular automaton [8|. 

In the present paper a finite version of Random 
Domino Automaton is investigated. The mathematical 
formulation in finite case is precise and the presented 
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results clarify which formulas are exact and allow to esti- 
mate approximations we impose in infinite case presented 
in [l|. We also show, that equations of finite RDA can 
reproduce results of |lj , when size N of the system is in- 
creasing and distributions satisfy an additional assump- 
tion (rii — > for big i). 

On the other hand, a time evolution of Finite RDA can 
exhibit a periodic-like behaviour (the assumption — >■ 
for big i is violated), which is a novel property. Thus, 
based on the same microscopic rules, depending on a 
choice of parameters of the model, a wide range of prop- 
erties is possible to obtain. In particular, such behaviour 
is interesting in the context of recurrence parameters of 
earthquakes (see e.g. |s l[Toj|). For other simple periodic- 
like models, see [111 

The finite case makes an opportunity to employ 
Markov chains techniques to analyse RDA. Investigat- 
ing the automaton in Markov chains framework we ar- 
rive at several novel conclusions, in particular related to 
expected waiting times for some specified behaviour. 

This article completes and substantially extends pre- 
vious studies of RDA on the level of mathematical struc- 
ture. We analyse properties of the automaton related 
to time evolution and others, as a preparation for fur- 
ther prospective comparisons with natural phenomena, 
including earthquakes. A matter of adjusting the model 
to the real data is left for the forthcoming paper. 

The plan of the article is as follows. Mimicking [l| in 
Section [II] we define the finite RDA. In Section Hill we de- 
rive respective equations for finite RDA. In Section ITVl wo 
will specify them for some chosen cases. In Section FVlwc 
will shortly describe Markov chains setting and describe 
time aspects of FRDA. Several examples are presented in 
Section [VT] The last Section fVIII contain conclusions and 
remarks. In the Appendix we show non existence of ex- 
act equations for RDA as well as present supplementary 
formulas and Table IXIVI displaying all states of RDA of 
size N = 10. 
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II. FINITE RDA 

The rules for Finite Random Domino Automaton are 
the same as in We assume: 

- space is 1-dimensional and discrete - consists of N cells; 

- periodic boundary conditions (the last cell is adjacent 
to the first one); 

- cell may be in one of two states: empty or occupied by 
a single ball; 

- time is discrete and in each time step an incoming ball 
hits one arbitrarily chosen cell (each cell is equally prob- 
able). 

The state of the automaton changes according to the 
following rule: 

• if the chosen cell is empty it becomes occupied with 
probability v\ with probability (1 — v) the incoming ball 
is rebounded and the state remains unchanged; 

• if the chosen cell is occupied, the incoming ball provokes 
an avalanche with probability p (it removes balls from hit 
cell and from all adjacent cells); with probability (1 — p) 
the incoming ball is rebounded and the state remains 
unchanged. 

The parameter v is assumed to be a constant but the 
parameter p is allowed to be a function of size of the hit 
cluster. The way in which the probability of removing a 
cluster depends on its size strongly influences evolution of 
the system and leads to various interesting properties, as 
presented in the following sections. We note in advance 
that in fact there is only one effective parameter p/v 
which affects properties of the automaton. Changing of p 
and v proportionally in a sense corresponds to a rescaling 
of time unit. 

A diagram shown below presents an automaton of size 
N = 12, with three clusters (of size 1, 2 and 4) in time t. 
An incoming ball provokes an relaxation of the size two, 
thus in time t + 1 there are two clusters (of size 1 and 4) . 



time = t 



time = t + 1 



The density p of the system is defined as 

JV 
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(3) 



In this article we investigate a stationary state of the 
automaton and hence the variables ni^nn,p and others 
are expected values and do not depend of time. 



III. EQUATIONS FOR FINITE RDA 

In this section we derive equations describing station- 
ary state of finite RDA. The general idea of the reasoning 
presented below is: the gain and loss terms balance one 
another. 



A. Balance of density p 

The density p may increase only if an empty cell be- 
comes occupied, and the gain per one time step is 1/N. 
It happens with probability ~ 2/(1— p). Density losses are 
realized by avalanches and may be of various size. The 
effective loss is a product of the size i of the avalanche 
and probability of its appearance pi(riii) / N . Any size i 
contribute, hence the balance of p reads 



1 N 

i=l 



(4) 



We emphasise, the above result is exact - no correlations 
were neglected. Its form is directly analogous to the re- 
spective formula in [l|. 



B. Balance of the total number of clusters 

Gain. A new cluster (can be of size 1 only) can be 
created in the interior of empty cluster of size > 3. 



Denote by n.;, i = 1, . . . , N the number of clusters of 
length i, and by n®, i — 1, . . . , N the number of empty 
clusters of length i. Due to periodic boundary conditions, 
the number of clusters is equal to the number of empty 
clusters in the lattice if two cases are excluded - when 
the lattice is full (single cluster of size N) and when the 
lattice is empty (single empty cluster of size N). Hence 
for 



N-l 



N-l 



n R =^2m, and n° R = ^ n- 



i=i 



we have 



(1) 



(2) 



(i— 2) cells = interior 



If the empty cluster is of size N, then each cell is in in- 
terior. Summing up contributions for all empty clusters, 
the probability is 



N-l 



i=3 



2\ n?i 



which can take a form (for N > 3) 



ft ^ i n R , "l 
V (l-p)-2u—+u—. 



(5) 



(0) 
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Loss. Two ways contribute: joining a cluster with 
another one and removing a cluster due to avalanche. 

Joining of two clusters can occur if there exists an 
empty cluster of length 1 between them. The exception 
is when the empty 1-cluster is the only one empty clus- 
ter, and the system consists of a single cluster of length 
N — 1. Hence, the probability of joining two clusters is 



The probability of avalanche is just 

N 

\ - riii 



i=l 



N 



(7) 



(8) 



Gathering these terms one obtains equation for balance 
of the total number of clusters n 



N(l -P)- 



N 

T.T, 



mi + n N -i = 2n R . 



(9) 



Again we emphasise that the above result is exact - no 
correlations were neglected. Finite size of the system 
reflects in the appearance of (2n R — njv-i) instead of In 
in the respective formula in |l| . 



C. Balance of n^s 

Loss. There are two modes, 
(a) Enlarging - an empty cluster on the edge of an i- 
cluster becomes occupied. There are two such empty 
clusters except for the case when system contains a single 
cluster of length N — 1. Hence, the respective rates are 



LV N 



N 



i = l,...,N 
i = N — 1. 



(10) 
(11) 



(b) Relaxation rate for any i — 1, . . . , N is given by 



(12) 



Gain. Again, there are two modes, 
(a) Enlarging. For N > 3, there are following rates de- 
pending on the size i of the cluster 



v(l-p)-2v!# + u%, 



2l/ — a T-l 



N 



1 = 1, (13) 

2 < i <N- 1, (14) 
i = N, (15) 



where a.E{i) is a probability that the size of empty cluster 
adjacent to the z-cluster is bigger than 1. It is clear that 



a N-2 



1 and 



0. 



(16) 



Formula (|T5j) does not have a factor 2, because there is 
only one empty cluster (of size 1). 



(b) Merger of two clusters up to the cluster of size i. 
Two clusters: one of size k G {1, 2, . . . , (i — 2)} and the 
other of size ((« — 1) — k) will be combined if the ball fills 
an empty cell between them. 




The probability is proportional to the number of empty 
1-clusters between fc-cluster and (i — 1 — fc)-cluster, 



"1 E 

P N 7i 



3<i<N-l, 



(17) 



where is a probability of such merger. For i — N there 
is a single cluster in the lattice (there are no two clusters 
to merge) - filling the gap between ends of (N — l)-cluster 
is already considered in (a). 

Gathering the terms, one obtains 



1 



7fl 



"2 



V 

1 

2^ + 2 

V 

1 



njv-i 



n N 



(N(l - p) - 2n R + n?) , (18) 
2niaf, (19) 
(2n i _ 1 af_ 1 +n? 7 f), (20) 

Miv-i^/ 1) | l ( 2 ^-2 + <7iY-l) , (21) 

BKN n N-U (22) 



where 3 < % < (N - 2). 

The last equation (l22l has simple explanation. The 
state with all cells being occupied (corresponding to n^) 
can be achieved only from the state with a single empty 
cell (corresponding to tin-i) with probability i/(l/N). 
On the other hand, the automaton leaves the state with 
all cells being occupied with probability pn- 

Note that equations (fI5]l and (|2"2"|) are exact. Correla- 
tions in the systems reflect in appearing of multipliers af 
and jf. Their values depends on possible configurations 
of states of the automaton. As shown in the Appendix, 
for TV > 5 exact formulas for af and as functions 
of riiS do not exist. Hence, it is necessary to propose 
approximated formulas. 



A mean field type approximation for af is 



E 



N-i 
k=l 



(23) 



For a given cluster of size i, the probability of appearance 
of an empty cluster of size 1 is calculated as proportional 
to the number of empty 1-clusters divided by the sum of 
the numbers of all empty clusters with size not exceeding 
N — 1, because there is no room for larger. 

When merger of two clusters up to a cluster of size i is 
considered, the room denoted by A is of size (N — 2 — (i — 
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l — k)) and the room denoted by B is of size (N — 2 — k) 
- see a diagram below. 



(i-l-fc) 




Hence a mean field type approximation for jf is of the 
form 



if « 7, A 



E ( v^JV- i'- 1- \-- ■ -i A-:--ji 



(24) 

It is also instructive to consider another approximation 



z-2 



7t ~7j 



EE 



fc=i 



n R n R 



(25) 



Section I VII contains quantitative estimation of proposed 
approximations. Comparison of this approximation with 
exact results for small sizes N is discussed in Section lVIII 



D. Thermodynamic limit 

In the paper [l[ an assumption of independence of clus- 
ters was considered. To have it adequate, it is required 
that there are no limitations in space, like those encoun- 
tered when formulas (l23l) and (j24l) were considered. For 
systems that are big enough, i.e., when N — > oo, an 
empty cluster adjacent to a given i-cluster can be of any 
size, and thus 



fe=l 



1 - -i 

n 



(26) 



This is consistent with the requirement that Hi — > 
when i — > oo, which is required to have moments of the 
riiS convergent. Similarly, 



i-2 



7(i) = E 



fe=i 



rik rii-i-k 



n 



ii 



(27) 



These formulas substituted into (|18l) - (j2"0)l give the respec- 
tive set of equations considered in The same reason- 
ing can be applied to balance equations. The form of 
equation Q is left unchanged under the limit. For equa- 
tion (JSJ), (2n R — njv_i) — > 2n, and it becomes of the 
form presented in [1]. 



IV. SPECIAL CASES 

For fixed form of rebound parameters equations de- 
scribing the automaton can be written in more specific 



form. This is the case for balance equations (j4| and ([9]), 
as well as for formulas for average cluster size 



/ •\ _ E£a ^ _ N P 



and average avalanche size 

V^JV -2 



(28) 



(29) 



We emphasize, these formulas are exact - correlations are 
encountered. We consider three special cases investigated 
in detail and illustrated by examples below. 



A. fi = const. 

For yU = const, and v = const, equation Q is of the 
form 



1/1^ -2 

na . 



and equation 



N{1- p(l + ^))+n N _ 1 =2n R . 



(30) 



(31) 



Also formulas for (i) and (w) are simplified only a little. 



B. — S/i where 9 — 5/u — const. 



Equation ((4]) is of the form 

(1 - p) = o P , 



(32) 



hence the density is given by remarkably neat (end exact) 
formula 



1 



P = 



1 



(33) 



Note that there is no dependence on the size of the system 
iV; for N — > oo it remains the same. 
Equation (0) can be written as 



N- 



1 



= (2 + 0)n R , 



(34) 



where we use equations (|2"2")l and (|3"3"|) . Hence the formula 
for n R is of the form 



n R = N- 



( u + i)(e + 2) (35) 

in direct analogy with n in N — > oo case Thus, 
n R plays the role of n, as indicated also in balance of n\ 
equation (fT8|) . The formula for n is 



n — n R + nj^ = N 



9(1 + e) n N 

77 ttt; — -r where e = . 

(0 + l)(6 + 2) n R 



(36) 
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The average cluster size is given by 

»-tt;(i+|)- <m 

The average avalanche size is equal to the average cluster 
size 

(w) = (i) , (38) 

because each cluster has the same probability to be re- 
moved from the lattice. 

The above formulas are exact (include correlations) 
and have good thermodynamic limit (e — > 0). Note 
also that variables p and ur depend on single parameter 
9. Formulas with dependence on 8 can be rewritten as 
functions of density p. 



state number 


example 11 


multiplicity 


contrib. to 


1 


^1 1 1 |«- 


1 




2 


->| | |.|<- 


3 




3 


^1 l-l-l^ 


3 


n 2 , n" 


4 


| • | • 1«| ^ 


1 


ri3 



Other states differ by translations. 



TABLE I: States for the size of the lattice N = 3. 



where assumed periodic boundary conditions are de- 
picted by hook-arrows. 

For periodic boundary conditions it is irrelevant to dis- 
tinguish between states which differ by a translation only. 
Hence, in example, we consider the following states equiv- 
alent: 



C. fj,(i) — r)/i 2 and \ — cr/v — const. 

Equation ^ is of the form 

N(l- p)= X {n R + n N ). 
Equation §§§ can be written as 



(39) 



1 N n 

N{l-p) = 2n R -x T n N +xY J - (40) 



where equation (|22|) is used, namely njv_i = Xjj n N- 
The average cluster size 



(i) = X-- 



1 ' 

and the average avalanche size 

1 



(w) = 



where 



(l + f)-(l--)f 



n N 



nn + n N e + 1 
Note that also these formulas are exact. 



(41) 



(42) 



(43) 



Thus states et^ are defined up to translational equivalence 
(see Tables HI and HI)) . The label numbers are assigned to 
the states, as shown in tables - no exact rule is applied. 

Further reduction of the number of states using reflec- 
tions can be done, but it is not very efficient procedure. 
We do not perform it, keeping symmetrical states sep- 
arate. They deliver a simple computation check - their 
probabilities are necessarily equal. 

Such space of states for the finite random domino au- 
tomaton is irreducible, aperiodic and recurrent. Transi- 
tion matrix P is defined by 



probability of transition 



For N = 3 the transition matrix is of the form 



(44) 



2/J2 

3 

V M3 



V 

Mi 
3 






'in 

3 





3 
3 





1 



(45) 



A*3 / 



where entries are found from analysis of transition prob- 
ability of all possible states a* (see Tab|T|. 
For N = 5 the transition matrix is 



V. FINITE RDA AS A MARKOV CHAIN 
A. General settings 

Finite Random Domino Automaton is a Markov chain, 
hence we use standard knowledge to solve several exam- 
ples for small N and derive a number of formulas for time 
aspects of the evolution of the system. 

In general, for the lattice of size N there are 2 N states, 



because each of N cells may be empty or occupied. 
N — 4, an exemplary state is 



For 



state number 


example* 


multiplicity 


contrib. to 


1 


■-►I ^ 


1 


ni 


2 


<-t | [ • | 


5 


ni,ni 


3 


t- H ' 1 • 1 • 1 ^ 


5 


n 2 ,n'i 


4 


^1 1 l-l l-l^ 


5 


m,ni,n2 


5 


<-»• I | • | • | • | 


5 


n 2 ,ri2 


6 


<— > | | • | | • | • | «-j 


5 


m,n2,ni 


7 


1 • 1 • 1 • 1 • 1 ^ 


5 


ri4 , ni 


8 




1 


n 5 



"Other states differ by translations. 

TABLE II: States for the size of the lattice N = 5. 
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pi 5 


- /xi - 4i/ 


2i/ 


2v 
















5 


- 2/Lt2 — 3^ 
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;/ 
















5 — 2/ti — 3v 




2;/ 








3/^3 











5 — 3m3 — 2v 





2v 








2/i2 


Mi 








5 - 2p 2 - Mi 


- 2v 2v 






















5 — 4/i 4 — v 


;/ 


V 5^5 




















5 - 5M5 



(46) 



Stationary distribution is given by 
v ■ P = v. 



(47) 



The number of states increase rapidly with N: for N = 
6 there are 14 states, for N = 7 there are 20 states and 
for N — 10 there are 108 states. The number of states 
for any N is bigger than 2 N /N, because translational 
symmetry of states is at most N, but always there are 
states with smaller symmetry, like empty state and fully 
occupied state. Thus practical usage of Markov chain 
settings for calculations is rather limited. This is one of 
the reasons for developing more "handy" framework, like 
presented in [l| and here. On the other hand, Markov 
chains can be used for illustrations and justifications of 
some properties, as presented below. 



B. Expected time of return 

As system evolves, it hits a given state many times. 
Here we consider expected value of the time of return 
from state with density p = to itself and next from the 
state with p = 1 to itself. 

Starting from state 1 (state with p = 0) the next state 
(different from state 1) contains a single 1-cluster only. 
This state - denoted by label 2 - has density p = 1/N. 
Expected time for this change is \ jv. 

Let Tj be the expected time to hit state 1 starting in 
state i. Then t\ = and for i =/= 1 

Ti = E(time to hit 1 | start in i) 

= l + J2pikE(l\k) = l + J2p ik T k , (48) 

k k 

where E(l|fc) = E(time to hit 1 | start in k). After solv- 
ing this system of equations, the return time is 



=1/v + r 2 . 



(49) 



Similarly, for state with p — 1 (state L) the next state 
(different from state L) is the empty state (with p = 0) 
and 



t L - 



i/mjv + h, 



(50) 



where fx is the expected time to hit state L starting in 
state 1. The respective equation to determine Ti for i =/= L 
reads 



PikTk, 



(51) 



and obviously tl = 0. 

Note that the expected time t^^L is equal to expected 
time of return from state 1 to state 1 through state L: 

th^L = h^L^x- (52) 
The expected time between two consecutive avalanches 



is 



(w) + 1 

1 - p. ' ' 



(53) 



where P r is the probability that the incoming ball is re- 
bounded both form empty or occupied cell: 



1 N 

P r = (1 - p)(l -v) + JjY. ^ l ~ I*)- 



(54) 



Note that (1 — P r ) is equal to the sum of probability of 
triggering an avalanche and probability that an empty 
cell becomes occupied, hence 



1 N 

P r + — ^ Mi™i* + (1 - P) v = 1- 



(55) 



i=X 



Formula ([531 can be derived as follows. In time between 
two consecutive avalanches, on average, (t av (l — P r ) — 1) 
cells become occupied in the system - it receives one ball 
per a time step, part of them are rebounded and one 
ball triggers the avalanche. An avalanche is reducing the 



probability of 


value 


rebound - occupied cell 
rebound - empty cell 
occupation of empty cell 
trigerring an avalanche 


££f=in*i(l-M0 
(l-p)(l-f) 

1 



TABLE III: Probabilities of all four possibilities occurring in 
a single time step during evolution of the automaton. 
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number of occupied cells by (w). These two quantities 
compensate each other, giving (|53l) . 

On the other hand, the expected time between two 
consecutive avalanches is equal to the inverse of the prob- 
ability of triggering an avalanche 



(56) 



Both expressions given in (|5"5|) and (|5"6"|) are equal to each 
other. 



C. Frequency distribution of avalanches 

The probability of states obtained from condition (|4"7} 
allows to determine the distribution of frequency of 
avalanches. The frequency fi of the avalanche of size 
i is given by the sum of products of probabilities of 
state k and respective transition probability pkj to the 
appropriate states j for all states that transition k — > j 
produce the avalanche of size i. 

For example, for N = 5, as can be seen in Table [TT1 
transitions 2 — > 1, 4 — > 2 and 6 — > 3 result in an 
avalanche of size 1, transitions 3 — > 1 and 6 — > 2 give 
an avalanche of size 2, transition 5 — > 1 of size 3, 7 — > 1 
of size 4 and 8 — !• 1 of size 5. Hence 



h 


= «2Ml/5 + 


i! 4 2/zi/5 + v 6 fj.i/5, 


(57) 


h 


= v 3 2^ 2 /5 - 


f u 6 2^2/5, 


(58) 


h 


= d 5 3^ 3 /5, 




(59) 


h 


= v 7 4^ 4 /5, 




(60) 


h 


= tW5, 




(61) 



where respective pkj are taken from transition matrix 
66]). 

The average time ti between two avalanches of size i is 
given by 



U = l/fi, 
in particular, for a maximum size N 

tL^L = In- 



(62) 



(63) 



The average time between (any) consecutive 
avalanches given by formula (|56p may be also cal- 
culated as 



(64) 



because the probability of avalanche of any size is just 
a sum of probabilities of all possible avalanches. In this 
way one can calculate also average time between any two 
consecutive avalanches of prescribed size - for example, 
size 4 and 5 (or any other subset of possible sizes). 




VI. EXAMPLES 

Below we present several examples to illustrate prop- 
erties of finite RDA as well as to demonstrate application 
of the schemes outlined above. 



N : 



This is the simplest non-trivial, worm-up example. For 
N ~ 3 the general results - i.e., for arbitrary Hi, [i 2 , /13 
and v - can be calculated explicitly. Usage of equations 
(fT8| - ((22|) leads to exact results as presented below (see 
Appendix) . The same can be also obtained from Markov 
chains framework. Equations (HJ, ^j, and (|2"2")l give 



where 



m 

n 2 = 3/D, 



D=— + ^i + 5^ + ^ + - 
2 2v v v A 



(65) 
(66) 
(67) 



M3 



From inspecting of Table U it is evident that n? = n 2 , 
n 2 = ni and n§ = 1 — n\ — n 2 — n 3 (all posibilities sum 
up to 1), hence 



n°=fl + ^ + 2^ + ^)/D. 
2v v v z 



General formulas for expected times of return are 

1 / 2^ 2 + 9fi 3 i> + 6fi 2 fi 3 
ti^i — — I 1 



v \ M3(Mi + 2v){2[i 2 + v) 

1 / v 11 jll 



tL^L = - [— + ^- + 5-+ 5 — 

v \[i 3 2 Iv v 



, 



The ratio Ai^i is 



2^2 



(68) 

(69) 
(70) 

(71) 



Note that it does not depend on fi 3 . If the probability of 
triggering an avalanche of size 1 and 2 is small compar- 
ing to the probability of occupation of an empty cell (i.e., 
/Ji/v w and [i 2 jv w 0) then ^l-s-l ~ ii->i- The next 
stage after the lattice is fully occupied is the empty state; 
hence, if these two average waiting times are comparable, 
then they occur with comparable frequency. That means 
quasi-periodic like behaviour of the system: within aver- 
age time 11 /2v the lattice become fully occupied, then 
the triggering of an avalanche of maximal size N occurs 
with average waiting time I/V3. The same can be ob- 
served for bigger sizes N. 
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-0.65 





,--« 1 


1 
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-0.54f 





FIG. 1: Plot of the Logio of riiS (left) and n°s (right) vs. i 
for N = 3 in three cases: pi = const, (dashed line), pi = 5/i 
(solid line) and pi = a/i 2 (dotted line). Rebound parameters 
are chosen to have density p = 1/2 in all cases (see main text 
for respective values). 



FIG. 3: Ratio of return times th-^L for N=5 (left) and 
N=7 (right) for three cases: pi — const, (top), pi — 5/i 
(middle) and pi = n/i 2 (bottom). Parameter t is equal to 
p/v, 5/v and n/v respectively. 



(i) 

H 



p, 



const. 



1.9281668 
2.2516538 



Uj = 5/i 



Pi = a/i 2 



2.1134407 
1.7226121 



TABLE IV: Average cluster size (i) and average avalanche size 
(w) for three different rebound parameters. Density p = 1/2, 
the size of the lattice N = 3. 



Figure [T] and Table HVl present examples of three types 
of dependence of rebound parameters on size i of clus- 
ters considered in Section IIV1 each having the same den- 
sity p = 1/2 (with 8 digits accuracy). To obtain this 
density we put for these three cases p/v = 0.444118 
(// = 0.444118, v = 1), e = 1 (6 = 1, v = 1) and 
X = 2.113440690 (77 = 1, v = 1/2.113440690) respec- 
tively. As seen from Figure Q] it is possible to obtain flat 
distribution for pi = 5/i - on that background, differ- 
ences between the cases are clearly visible: pi = const. 
discriminate the existence of big clusters fostering big 
avalanches; the opposite is for pi = a/i 2 . Average clus- 
ter size and avalanche size data presented in Table IIVI 
confirm this conclusion. 





FIG. 2: Plot of the Logio of nts (left) and n° s (right) versus i 
for N = 5 in three cases: pi = const, (dashed line), pi = 5/i 
(solid line) and pi = a/i 2 (dotted line). Rebound parameters 
are chosen to have density p = 1/4 in all cases (see main text 
for respective values). 





Pi — const. 


Pi — 5/i 


Pi = a/i 2 


<*> 


1.427017126 


1.632218845 


1.985611461 


(w) 


1.845355789 


1.632218845 


1.41360643 



TABLE V: Average cluster size (i) and average avalanche size 
(w) for three different rebound parameters. Density p = 1/4, 
the size of the lattice N = 5. 



B. N = 5 

For TV = 5 it is impossible to write down exact equa- 
tions (fl~g)) - (f2"2l depending on values of 71^ only - see 
Appendix for details. The case can be solved as a 
Markov process, but obtained general formulas are rel- 
atively complicated. 

In this example we investigate properties of the system 
with density p = 1/4. Figure [2] and Table fVl compare 
results in three cases: p^/v = 16257/10000 the density 
p = 0.2500003184; for 9 = 3 the density p = 0.25 exactly; 
and x = 5.95682 gives the density p = 0.2500004527. 

General expressions for return times fi_».i and as 
well as their ratio (presented in Appendix) are relatively 
complex. Note that the return times - except of the de- 
pendence on t - are proportional to \/v. Below we spec- 
ify the ratio iz,->-L Ai-s-i in three cases: for pi = const, 
where t = pi/v, it is equal to 



24i fa 



154i 5 + 413i 4 + 586^ + 467i 2 + 182i + 24 



, (72) 



24£ 2 + 54t + 24 
for pi — 8/i, where 6 = const and t = S/v, it is equal to 

At 6 + 40i 5 + 169t 4 + 395£ 3 + 550t 2 + 432i + 144 

56i 2 + 192i + 144 ' ^ ' 

and for pi = a/i 2 , where a — const and t = a /v, is 

2t 6 + 39i 5 + 304i 4 + 1232i 3 + 2840i 2 + 3744i + 2304 
496i 2 + 2208t + 2304 ' 

In each case the ratio is a rational function of t, which 
is equal to 1 for t = and asymptotically ~ i 4 for t — > 
00. A generalisation of this observation is a Conjecture 
formulated in Section I VIII A comparison of these ratios 
is presented in left part of Figure |3l Table I VII shows 
that for the cases discussed above with average density 





Pi = const. 


Pi = 5/i 


Pi = a/i 2 


R = tL~+L /tl^l 


w 52.212 


m 35.441 


w 34.801 



TABLE VI: Coefficient R = th^h/tx^x for three different 
rebound parameters (see main text for details). Density for 
all cases p = 1/4, the size of the lattice N = 5. 



123451 2345 



1234567 1234567 



FIG. 4: Left. Plot of Logio of Us vs. i for three rebound 
parameters for fixed density p = 1/4 for TV — 5. Three cases: 
/ii = const, (dashed line), fj,i — 5/i (solid line) and pi — a/i 2 
(dotted line). Rebound parameters are chosen to have den- 
sity p — 1/4 in all cases (see main text for respective values). 
Right. Plot of Logio of ijS versus i for various densities for 
rebound parameter of the form pi — 5/i for TV = 5. Densi- 
ties are chosen as rL, j, |, |, thinner line corresponds to 
smaller density. 

/0=l/4 the highest value of i? is for /x^ = const, and the 
smallest for /ij = a/i 2 (not much different from the value 
for /ii = S/i). 

Average waiting times U for avalanche of size i can be 
also found. For example for /i; = S/i, where 5 — const., 
they are presented in the Appendix (equations (|81[1 - 
(1851) ). The average time between any two consecutive 
avalanches is 

_ 4t 5 + 48t 4 + 237t 3 + 603i 2 + 762t + 360 
tav ~ z4(4i 4 + 36i 3 + 121t 2 + 168i + 72) ' ( ' 

where t = 5/ v. All these quantities arc proportional to 
1/za Figure [4] in the left panel presents waiting times 
ti in for fixed density p = 1/4 in three cases mentioned 
above. There are no big differences both in character of 
dependence of ti on i and also values of t av do not differ 
much: for Ui = const, average time is t av 24.60, for 
Hi = 5/i it is 21.76 and for pn — a/i 2 it is « 18.85. 
(Choosing parameters to have density p = 1/4 we put 
v = 1/10 for all cases.) 

Average waiting times ti, i — 1, . . . , 5 in the case pi — 
S/i for various densities are shown in the right panel of 
Figure |U For small densities the maximal waiting time 
ti is for i — 5, while for bigger densities the maximum 
is for i = 3. Average waiting times range from k, 13.57 
for p = 1/10 through « 21.76, « 50.22, w 145.01 for 
densities 1/4, 1/2, 3/4 respectively, up to ~ 441.60 for 
density p = 9/10. (Again v = 1/10 for all cases.) 

C. N — 7 

For N = 7 we investigate properties of the system 
with the density p = 3/4. Parameters are chosen as 
follows: fi = 1, v = 173024/10000 gives the density 
p = 0.7500001621, 9 = 1/3 gives p = 3/4 exactly, and 
(i = 1, v := 1000000/1578886 gives p = 0.7500002817. 
Distributions of clusters are presented in Figure [S] and 
average cluster and avalanche sizes in Table IVIII Again 
differences in distributions are not big, but average 



FIG. 5: Plot of the Logw of rus (left) and n°s (right) versus i 
for TV = 7 in three cases: pt = const, (dashed line), pi = 5/i 
(solid line) and pt = a/i 2 (dotted line). Rebound parameters 
are chosen to have density p = 3/4 in all cases (see main text 
for respective values). 

avalanche size differs significantly between considered 
cases. 

The novel property visible in the figure is that the high- 
est probability is for the cluster of maximal size i = TV . 
Thus, the system prefers merging clusters for high den- 
sity. 

A comparison of the ratios of return times R = 
tL^L/ti^i is presented in the right panel of Figure [31 
while formulas are presented in the Appendix. In each 
case the ratio is a rational function of t, which is equal to 
1 for t = and asymptotically ~ t e for t — > oo, which 
supports a Conjecture formulated in Section IVlTl Table 
IVIIII shows that for the cases discussed above, with av- 
erage density p = 3/4, the highest value of the ratio R is 
for p,i = 5/i and the smallest for pi = const, (which does 
not differ much from the value for p,i = 5/i). This is an 
opposite order comparing to the case with p = 1/5 for 
TV = 5 considered above. Thus, for higher densities the 
automaton prefers more periodic-like behaviour when it 
is relatively easier to trigger big avalanches. 

The size TV = 7 is big enough to notice how 
the actual density of the system (possible values are 
Q,ip,f,f,f,f,f,l)is distributed for various average den- 
sities. Results are shown in Figure |51 For small densities, 
like p = 0.2, the maximum is for small i, that means that 
big densities and big avalanches are rare. Then, when the 
density increases, the bell-like shape distribution appears 
and its maximum is shifted to the bigger values. Next, 
for densities like 0.6 or bigger, the maximum probability 
is for biggest possible size i = TV and the most probable 
state is that with p = 1. To achieve big average density, 
the system must spend a substantial time being fully oc- 
cupied. The evolution of such a system consists of two 
phases: filing up and waiting for avalanche of maximal 
size, as is described above while discussing the times of 
return for TV = 3. 





Pi = const. 


Pi = 5/i 


Pi = a/i 


<*> 


4.274328495 


4.385371765 


4.736665115 




5.767461682 


4.385371765 


2.671314107 



TABLE VII: Average cluster size (i) and average avalanche 
size (w) for three different rebound parameters. Density p = 
3/4, the size of the lattice TV = 7. 
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FIG. 6: Probability distributions of actual density of the 
system for the case (Xi/ V — 9/i for various average densi- 
ties of 0.2, 0.4, 0.5, 0.6, 0.8 - the respective parameters are 



9 = 4, ■ 



2 1 
3' 4- 



On the plot: smaller average density cor- 



responds to the higher probability rate for density equal to 
(and to the lower probability rate for density equal to 1). 
The size of the system is N = 7. 



For N — 500 and constant parameters fj, = v = 1, nu- 
merical experiments show that the density fits a Gaussian 
distribution [13[ . 



D. 



N = 10 



In the example with the biggest N presented here we 
investigate in several cases influence of correlations and 
compare exact results with proposed approximations for 
af, and 7^*. On the other hand, size N = 10 
requires relatively complex calculations - the transition 
matrix is of size 108 x 108 and has about 1000 non-zero 
entries. All possible states are presented in Table IXJVl 

The size N — 10 is the smallest with states which con- 
sist of the same clusters, but in essentially different order. 
(For smaller N states with different order of clusters were 
equivalent with respect to reflections.) Namely, the state 



and the state 89 



In this subsection we consider also the relative differ- 
ence between probabilities of these two states, namely 
A = (.P88 — 2p89)/P88 for various rebound parameters as 





[ii = const. 


Hi = S/i 


fli = a/i 2 


R = ti-fi/tl-H 


« 1.4844 


» 1.6887 


« 2.7001 



TABLE VIII: Coefficient R = th^h/t\^\ for three different 
rebound parameters (see main text for details). Density for 
all cases p = 3/4, the size of the lattice N = 7. 
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FIG. 7: Plot of the Log 1Q of n^s (left) and n° (right) for 
N = 10 in three cases: fii = const, (dashed line), Ui — S/i 
(solid line) and /it — a/i 2 (dotted line). 



a measure of adequacy of independence of clusters as- 
sumption. The multiplier 2 in the above formula is nec- 
essary because the multiplicity of state 89 is equal to five, 
and the multiplicity of the state 88 is equal to ten. This 
quantity reflects dependence of respective probabilities 
on specific order of clusters in the system. We assume 
there is no such dependence in order to write down ap- 
proximations af, jf and "if R . 

Other quantities analysed in examples below are 

n E - n A n E -n A -v E --y A -y E -~, AR -y E -~y A , -v E -■y AR 

a l a l a 4 a 4 75 ~5 15 75 79 79 arid 79 7g 



af ' af 1 ' 7e 7g 7s 

These quantities measure the quality of approximation 
formulas for i = 1 and 4 for a coefficients, and for i = 5 
and 9 for 7 coefficients - just to test approximations for 
n.2 — the first approximate equation, n$ - the middle one, 
and rig the last one (those for ri\ and n w are exact). 
A formula for exact value of af , obtained from detailed 
analysis of states of the automaton, is presented in Ap- 
pendix. 

a. Cases with constants equal to I. As a first set we 
consider three cases with the minimal possible rebounds 
factors, i.e., we put all constants equal to 1. Cases with 
/Xj = 1, /Xj = 1/i and ^ = 1/i 2 with v = 1 are presented 
in Figure and Table ED 

Three different rebound parameter types result in var- 
ious average density values, and hence different distribu- 



P 

<*•> 
(w) 

PSS ~2p 89 



7 5 -7 5 







AR 

JS 

~E 



7 5 
79 E -79^ 



79^ 



m = 1 


m = i/i 


Hi = 1/i 2 


0.3076370614 


0.5 


0.8822697788 


1.5985438 


2.872872532 


7.493995763 


2.250583644 


2.872872532 


3.725820785 


0.00865 


0.01899 


0.01868 


0.00909 


0.066784 


0.19468 


0.08795 


0.07428 


0.09239 


-0.01141 


-0.28676 


-0.89662 


-0.00110 


-0.06842 


-0.13292 


0.35717 


0.29878 


0.17045 


0.48970 


0.62573 


0.71428 


18.51 


25.59 


96.28 


379.61 


5.2988 


1.4606 



79 



TABLE IX: Three cases: a = 1, 8 = 1 and a 
v = 1) for the size of the lattice iV = 10. 
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FIG. 8: Plot of the Logio of n^s and vs. Log\o{i) for 
N = 10 in two cases: pi = const, (solid line) and pi = S/i 
(dashed line) - upper panels, and in case pt = 5/i 2 - lower 
panels. 



tions. In all cases, the assumption of independence of 
clusters is well satisfied; the respective error A does not 
exceed 2%. An approximation for af is less than 10% 
for all cases, but af strongly depends on the case (in 
fact it depends on density, as will be seen below). Ap- 
proximation formulas for je perform in diversified way - 
jar is better for mid i terms, while ja is better for big 
i terms. Nevertheless, both cases provide rather roughly 
appropriate values. These examples also suggest that for 
higher densities the system exhibits a periodic-like evo- 
lution. 

b. Big densities In order to investigate evolution of 
the system with high average density (and strong devia- 
tions in actual density) we consider case p = const with 
pi = 1/100 and v = 1, which gives the density p k, 0.91, 





m ~ 1 


m ~ i/j 


Pi ~ 1/i 2 


p 


0.9145269069 


0.9145269069 


0.9897960692 


(<) 


7.764600567 


7.805017612 


9.700144892 


(w) 


9.346154002 


7.805017612 


8.33021261 


P88— 2p 89 
£™ A 


0.00161 
0.30898 


0.00464 
0.28421 


0.00191 
0.31377 




-0.01511 


0.01520 


0.00330 


7,?-7jf 
1 E 


-0.96763 


-0.93825 


-1.05131 


-0.15298 


-0.14354 


-0.14789 


7f 


0.42749 
0.80211 


0.39218 
0.78888 


0.38965 
0.80019 


119.18 


123.60° 


1002.44 


th->L /t\^\ 


1.1066 


1.1339 


1.0277 



"The system stays in fully occupied state l//xio ~ 107,5, which 
is longer than 100 as in fi = const, case (previous column). The 
respective average times for filling up the lattice are 19 and m 16. 



TABLE X: Three cases with "big" p for the size of the lattice 
N = 10 (see main text for details). 



and case S/i with p x = 4673077001/5 * 10 10 « 0.093 and 
v = 1 to obtain the same density (with 10 digits accu- 
racy) for comparison. Also we consider case of a/i 2 with 
pi = 1/10 and v = 1 which gives the density p ss 0.99. 
The results are presented in Figure [8] and Table [X] 

Plots of respective distributions for p t ~ 1/i and pi ~ 
1 are overlapping each other. For relatively small size 
N = 10, fixing the average density of the system strongly 
determines distributions, making the dependence on re- 
bound parameters not essential. Their influence becomes 
more visible for larger sizes N of the lattice. In case of 
high density, the system just spend much time being fully 
occupied. 

For high densities, the assumption of independence of 
clusters is well satisfied; the respective error A does not 
exceed 0.5%. An approximation for af is fairly good 
(w 1.5% or less), but af has only accuracy ~ 30%. Ap- 
proximation formula for jar is much better for mid i 
terms (though giving only « 15% accuracy), while ja is 
better for big i terms (« 40%). Thus, for high density 
cases the proposed set of equations for n^s does not re- 
produce actual distribution. Note, however, that there 
are other exact equations valid for any density. 

The parameter Ai^i for /i^ ~ 1/i case is bigger 
than for pi ~ l case (both cases have the same "big" 
density), which agrees with the results for N — 7 with 
p = 3/4 presented in Table EEID 

c. Small densities To present system behaviour in 
small average density we choose p± = 1 and v = 1/10 for 
case p — const - it gives density p w 0.08. Then for the 
remaining two cases we have the same density p w 0.01 
(with 10 digits accuracy), with the following parameters: 
pi = 1 and v = 50000000/4798952601 w 0.01 - for case 
S/i and p\ = 1 and v = 1/100 for case a/i 2 . The results 
are presented in Figure [9] and Table IXII 

For small densities assumption of independence of clus- 
ters is well satisfied. In general, all proposed approxima- 





Pi ~ 1 


Pi ~ 1/i 


Pi ~ 1/i 2 


p 


0.0779280356 


0.01031150521 


0.01031150521 


<i> 


1.09149321 


1.02083788 


1.041894016 


(w) 


1.183235221 


1.02083788 


1.020408163 


P88— 2p s9 
I 88 A 


6.9231 * 10~ 5 
1.7483 * 10" 5 


6.4609 * 10~ 4 
2.7693 * 10" 7 


3.5663 * 10^ 3 
1.0585 * 10~ 6 


E A 


0.14229 


0.16313 


0.16274 


7.f-7 5 A 
75 K 

7 5 B 


0.0009069 
0.0009101 


0.0020519158 
0.0020519262 


0.0040100 
0.0040124 


7 9 E -7# 
7<f 

7 9 B 


0.09615 
0.12138 


0.01473 
0.02356 


0.00606 
0.03305 


22.32 


106.35 


110.57 


tL->L/tl->l 


2220903488.0 


1.666292752 • 10 14 


9971770329.0 


~ 2 ■ 10 10 


~ 2 ■ 10 14 


~ 1 • 10 11 



TABLE XI: Three cases with "small" p for the size of the 
lattice iV = 10 (see main text for details). 
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2468 10 2468 10 



FIG. 9: Plot of the Logio of n^s and n° for TV = 10 in case 
Pi = const. - upper line, and in two cases: pi = S/i (solid 
line) and pi = S/i 2 (dashed line) - lower line. 

tions are fairly good. An approximation for af is the 
worst; its accuracy is only w 15%. As previously, ap- 
proximation formula for 7^4 is better than ^ar for big 
i terms, but it appears that for mid i terms both for- 
mulas give almost the same values (because n^s decrease 
rapidly). Thus, for small densities the set of equations 
for riiS can be used to reproduce the actual distribution. 

It is very improbable to find the lattice fully occupied 
for small average densities, which is reflected in high val- 
ues of the parameter R = £l->lAi->i- The parameter 
R for fii ~ 1/i case is bigger than its for pi ~ 1/i 2 case 
(both cases have the same " small" density) , which agrees 
with the results for N — 5 with p = 1/4 presented in Ta- 
bleEH 



VII. CONCLUSIONS 

In this article we investigated in detail a finite ver- 
sion of one-dimensional non-equilibrium dynamical sys- 
tem - Random Domino Automaton. It is a simple, slowly 
driven system with avalanches. The advantage of RDA 
(comparing to Drossel-Schwabl model) is the dependence 
of rebound parameters on the size of a cluster. This 
crucial extension allows for producing a wider class of 
distributions by the automaton, as well as leads to sev- 
eral exact formulas. Exponential type and inverse-power 
type distributions of clusters were studied in [l|; the 
present work examines also V-shapc distributions and 
quasi-periodic like behavior. 

Detailed analysis of finite RDA, including finite size 
effects, extends and explains the previously obtained re- 
sults for RDA. Moreover, we also analyzed approxima- 
tions made when deriving equations for the stationary 
state of the automaton. This allows for the following 
conclusions. 

The balance of p equation (0| and the balance of N 
equations (jH]) are exact - their forms incorporate all cor- 
relations present in the system. The first one has a form 



independent of the size of the lattice N, thus it is exactly 
the same as for RDA. The second one contains correc- 
tion for finite size effect, namely a term (2riR — njy-i), 
which replaces the term 2n for RDA. When 7ijv-i and 
tin are negligible, these two terms coincide. For finite 
RDA, balance of n^s equations (fT5)) - (|!?!2)) contains two ex- 
tra equatins, for i — N — 1 and i = N, comparing to the 
those for RDA. The first (for ni) and the last (for rijv) 
are exact. Note that all those equations are written for 
rebound parameter p = p(i) being a function of cluster 
size and v being a constant. 

The most remarkable special case is when p = S/i, 
when any cluster has the same probability to be removed 
as an avalanche independently of its size i. It appears 
that the system depends on a single parameter 9 = 5 /v, 
or equivalently, due to neat exact formula (eq. (|3"3")l ) 

p = TTe> 

the properties of the system may be characterized by 
the value of the average density. Note that the above 
expression does not depends on the size N, and is the 
same as for RDA. This specialization leads to more neat 
formulas, like the equation for tir (eq. (|35[) ) 



Note again that it has the same form as for RDA, except 
that n is replaced by riR (n = ur + n^). Summarizing, 
the model allows to derive a number of explicit depen- 
dencies, as shown in Sections IIIII and HVl 

The Random Domino Automaton defines a discrete 
time Markov process of order 1 and, in principle, may 
be solved exactly. However, it turns out that computa- 
tions are fairly complex and exact formulas are long, as 
visible from examples presented in the Appendix. Also, 
the exact numerical values are in the form of big num- 
bers — in every considered example (N — 3, 4, 5, 6, 7, 10) 
significantly big prime numbers were encountered. For 
example, for the simplest possible rebound parameters 
(p = 1 and v = 1) the exact value of denominators of 
probabilities of states for N = 10 (see Table IXIV[) is 
a 65-digit integer. Its prime factorization (presented in 
the Appendix) contains a 56-digit integer, which cannot 
be simplified with numerators. Thus, the usefulness of 
Markov chains for finding both formulas and values n^s 
is limited in practice. 

Nevertheless, Markov chains framework leads to inter- 
esting results concerning analysis of times of recurrence 
for specific states. A return time to the state with density 
p = 1 (equation ([50)1 ) 

tL->L = + n, 

consists of two parts: waiting of fully occupied lattice for 
triggering a maximal avalanche and "loading" time, when 
the lattice is filled up, respectively. If the average density 
of the system is small, the second time is very long. The 
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formula is more interesting for systems with relatively big 
average density, when the "loading" time is comparable 
to waiting time for triggering the biggest avalanche. Such 
a system exhibits a periodic like behavior. Dividing the 
waiting time by the waiting time ti_>i (given by 

equation (|49l) ) one has the following measure of quasi- 
periodicity 

R = tL->L/tl->.l = il-s-L-s-lAl-s-l- 

If R = 1 then the system is periodic. 

Several considered examples lead to the following con- 
jecture concerning the coefficient R. 

Conjecture. The ratio of return times s-lAi >i as 

a function of t being the ratio of constants from rebound 
parameters (fi/v, S/v, <j/v) are rational functions of t, 

f(t) = tjj s-l/^i s-i with the following properties for 

any size A of the system 

f(t = 0) = 1, (76) 
f(t) 

lim = const. (77) 

t— s-oo t 1 

The conjecture relates the size of the system A with 
asymptotic behavior of ratio of waiting times. 

There are big fluctuations (variations of actual den- 
sity) during the evolution of systems with relatively big 
average densities. If the system is likely to achieve a fully 
occupied state, the next state is an empty state, and the 
variations in density are maximal. Nevertheless, some 
parameters of stationary state (more precisely, statisti- 
cally stationary state) satisfy exact equations, as shown 
above. For big average densities, the system fluctuates 
within the whole possible range, and cannot be thought 
of as having approximately stationary values durin g th e 
evolution. This aspect is easy to be overlooked (see 

It is argued in the Appendix that no exact equations 
for riiS exist for the size A > 5. Thus, to have compact 
equations for ms, some approximation formulas are pro- 
posed. The first general conclusion from the examples is 
that the approximations are acceptable for small densi- 
ties, but for big densities the errors are substantial. The 
main reason is that for big densities correlations become 
more important and fluctuations makes actual values of 
the parameters substantially different from their station- 
ary values, which are present in the formulas. These 
properties are particularly severe for small sizes of the 
system, where every avalanche changes the actual den- 
sity considerably. 

Table IXIII presents a dependence of a relative error of 
af with respect to af on size N of the system. For big- 
ger A the accuracy of approximation is growing, which 
corresponds well with the remark in the last paragraph. 

It can be noticed from the distributions of n^s of ex- 
amples presented above that all n^s except of the last 
two (namely Un-i and npf) are placed on one "regular" 
curve, while the last two deviate from it. It may be re- 
garded as a (correction of) finite size effect. Also in the 
respective set of equations (|18 p -(|22 p . the last two (for 
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1.33(3) 
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« 0.6316 


w 0.053 


0.32(32) 
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~ U.DODO 
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r^i u.oioy 


6 


« 0.685296 


w 0.669232 


w 0.023 


w 0.3102 
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« 0.685523 


« 0.675066 


w 0.015 


w 0.3086 


10 


« 0.685436 


w 0.679205" 


w 0.0091 


w 0.3076 


4000 6 


« 0.677 


« 0.677 




« 0.3076 



"(1 -nj/n) 0.679229 
6 Results from Q 

TABLE XII: Comparision of af with for N = 3, 4, 5, 6, 7 
and 10 for parameters fi = 1 and i/ = 1. The last line presents 
value of af for iV = 4000 obtained from simulations and 
equations (respectively) in 0|. 

i — A — 1 and i = A) have a form different from the pre- 
vious ones. Thus, neglecting the size restriction, which in 
fact ignores the last two equations, is justified when the 
deviations of the last two n^s from the "regular" curve 
are not big. That happens for small densities. 

It appears also that for index i in his middle range of 
values an approximation formula jar works better than 
7>1, in spite of the fact that it looks to be more rough 
approximation. For distribution of rijS vanishing rapidly 
(i.e., for small densities) both give comparable results. 

All this justifies the form of equations for n^s presented 
in [l| as valid for small densities. A detailed examination 
of the RDA for big densities requires further investiga- 
tions. 

This article explores properties of FRDA in order pre- 
pare to modeling of real data. In this context, among 
others, formulas for waiting times can be used. We epha- 
size also a formula ([55]l 

_ (w) + l 
av ~ 1 - P r ' 

which relates the measure of scattering (dissipation) of 
balls P r with the average size of avalanche (w) and av- 
erage time between any two consecutive avalanches t av , 
which are a priori measurable quantities. 

The Random Domino Automaton proved to be a 
stochastic dynamical system with interesting mathemat- 
ical structure. It may be viewed as extension of Drossel- 
Schwabl model, and we showed that this is a substantial 
generalization with a wide range of novel properties. W 
expect it can also be applied to natural phenomena, in- 
cluding earthquakes and forest-fires. This is our aim for 
the future work. 



Appendix. Exact equations for iV = 3, 4 and their 
non-existence for n > 5 

For arbitrary size A, there are four exact equations: 
balance of p - equation @ , balance of n - equation , 
for m - equation ([18]) and for - equation ([22]). 
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state number 
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multiplicity 


contrib. to 
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^ 1 1 1 1 \ 4r^> 
^ 1 1 1 ^ 




n Y 


2 


C_ ^ 1 * ^ 


4 


'"1 ) '"J, 


3 


^ 1 1 * 1 * ^ 


4 


n 3 ,ri2 


4 


| | • | | • ^ 


2 


ni , n" 


5 


^ • | • | • | <-J 


4 


n 3 ,ni J 


6 


■— > | • | • | • | • | ^ 


1 


ri4 



"Other states differ by shifts. 

TABLE XIII: States for the size of the lattice N = 4. 

Size N — 3. Equation for ri2 is of the form (|2ip. namely 
1 



"■2 



2^ + 1 



(2maf + n? 72 B ) 



In this case the only companion to single one-cluster is 
an empty two-cluster (see state 2 in Tab|I|, hence 

af = 1 and 7 f = 0. 

Thus, we arrived at the exact form of the equation for 
n 2 . 

Size TV — 4. All states of the automaton and their labels 
are presented in Tab lXIlil Equation for n 2 is of the form 
dH 



V 

where to a E contributes only state 2, not state 4. Hence 

2P4 



P2 
P2 + 2p4 



= 11- 



= 11- 



rVj - n 3 
m 



where pi is probability of state i. Thus af is expressed 
as function of n^s and n?. The equation for is of the 
form (121]) 



1 



nz = 



- (2n 2 a» + n^i) 



In this case 



of = 1, 



J 



because only state 3 contributes. The state 4 (and not 
state 5) contributes to jf , therefore 



73* 



2 Pi +p 5 I n°J \ n\ 



This completes the task of writing exact equations for 
N = 4. 

Size N = 5. States and their labels are presented in 
TabHH In this case, the coefficients are as follows 



af 



P2 + Pi 

m 







P3 E Pa 

— ) 74 = — 
n 2 



n 



Pi - " 4 



n 







Summing up the probabilities contributing to n\,ni and 
n 2 one obtains 



n? = _p 4 

"1 = P2 

n 2 = P3 



-2p 6 
-2 Pi 
Pg- 



■ «4, 
P6, 



(78) 
(79) 
(80) 



The set cannot be solved for p 2 ,P3,Pi,P6- Since there 
are no more equations for those coefficients, respective 
as and 7s cannot be expressed as functions of riiS only 
in an exact manner. 

Sizes bigger than 5. An argument for non-existence 
of exact set of equations (jT5]) - (|2"2"|) . i.e., non-existence of 
exact formulas for af s and 7^s as functions of rijS and 
n\ is based on the same impossibility of solving equations 
as presented above. 

An increase of size of a grid N by 1 results in an in- 
crease of the set of n\,n-z,. .. by one and much bigger 
increase of the number of states. An analog of the set 
of equations (|78|) - (f5U|) will contain much more probabil- 
ities of states pi,p 2 , ■ ■ ■, on the right hand side - there 
will be more states containing 1-clusters, 2-clusters and 
so on, and contributing to ni,n 2 , ■ ■ ■ respectively. Thus, 
it is impossible to express those probabilities of states as 
functions of n\, m, n 2l . . . only. As a consequence, there 
are no general exact formulas for afs, and jf. 



Appendix. Formulas 

The return times for N = 5 for general values of the parameters /ii, fi 2 , p,$, /Z4, /Z5 and v are 



(48/^ 5 + 
560/ii/X5^ 4 



v ' /j 5 (4/j 4 +;y)(3/J3+2i>)(2mu : 1 /J2+2M?i'+4M^2 + 17 ' ^i^2V+l3^{v 2 + li^!vlv+i3^i^2V 2 +33^iv 3 + l&filv 2 +i2fi 2 v 3 +3&V 4 -) 



60/i 2 ^ 5 + 54ti 3 ^ 5 + 750^ 5 + 72z/ e + 8^ + 8^ 4 + 20^iM2^ 4 + 12/n^ 4 + 24^3^ - 
+- 760ri 2 /i 5 ^ 4 + 720^3^5^ 4 + 1560^4A*5^ 4 + 100^/i 5 z^ 3 -I- 160/i 2 .^ 3 + AQ^l^v 2 + 2Q^ 2 ^v 2 - 
60 fil[i 3 [i 5 is 2 + 240 nlfiifisis 2 + 150/i2^3A*5^ 2 + 480^2^4^5 1,2 + 340^iii 2 ^5^ 3 + 390^i^3^ 5 ^ 3 + 1280^i/i4^ 5 ^ 3 - 

645^ 2 ^ 3 it 5 ^ 3 + 1840^2^4^ 5 ^ 3 + 1800^ 3 ^4/Z 5 ^ 3 + 240iii^2M3M4A t 5 + 120^1 2 y U 2 A*3A t 4M5 + 285^!tt2^3^ 5 ^ 2 +60^i/^| ^3^5^" 
30/i 2 /i2/^3M5^ +960/ii/i2/i4/^5^ 2 + 160/ii/i2M4M5^ + 80/i 2 /i2/X4/Lt5l'-|-1320/ii/i3/i4 
600^3^4/^ + H40/ii^2/Lt3/X4M5^), 
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th^L = — + o WU ixm UTa TJj ,1 2 , on TTo IS 2.07 T x (626/ii^ 5 + 844u 2 ^ 5 + 828u 3 ^ 5 + 1848U4^ 5 + 

822v 6 + 126/zfz/ 1 + 4/^ 3 + 192^ 4 + 8/i 2 ^ 2 .^ 2 + 426^i^ 4 + 489/ii/i 3 f 4 + 1544^ 4 ^ 4 + 771/i 2 M3^ 4 + 2176^ 2 M4^ 4 + 
2232^ 3 ^ 4 +68/xi/z!i/ 3 +54/Lt^2^ 3 +4/i?^2^ 2 ^ 

48/i 2 /i 2 M3M4 + 102^i^2M3^ 2 + 81miM2A*3 1/2 + 12m 2 M2M3^ + 272u\u\iiav 2 + 216^lfX2^i^ 2 + 32^l^U4iy + 396^/i 3 /i4^ 2 + 
792^2M3/M^ 2 + 24/xf;«2M3M4 + 414^i/x 2 M3i /3 + &n\^2^v + 1304^i^ 2 /"4^ 3 + 16/if/xa/X4^ + 1716^i/i 3 ^4^ 3 + 2Au\azH 4 v + 

2604/i2Ai3M4^ 3 + 1656^i//2M3M4^ 2 + 408^i^2M3M4^ + 324/1^2^3^4^) ■ 

Their ratio is 

, _ (4^4+i/)*(3M3+2t/)*(2/jf/j 2 +2M? ^+4/j| ^2 + 17^1^2 t/+13/jf i^ + M/jip 2 , 1^+43^1 M2^ 2 +33/ji 1^ + 16^2 ^ 2 +42ai2^ 3 +36<>) 

EL^L/tl^l — 2^(4 Al 5 ! + 10piM2+24Mii'+6M3A'i+4M2+30p2^+12M3M2+36i/ 2 +27A 1 3i/) ' 

Average waiting times ^ for avalanche of size i in case /ii = <5/i, where (5 = const, and £ := <5/z/, for N = 5, are 

4i 5 + 48i 4 + 237i 3 + 603i 2 + 762i + 360 
<x ~ i4 2 (4i 3 + 28i 2 + 69i + 60) ' ( 3 



4t 5 + 48i 4 + 237t 3 + 603i 2 + 762< + 360 
2^ 2 (4i 2 + 16t + 15) : 



(82) 



At 6 + 56i 5 + 333i 4 + 1077i 3 + 1968t 2 + 1884i + 720 
h ~ 2^ 2 (10i 2 + 31i+18) ' ( ' 

At 7 + 60£ 6 + 389t 5 + 1410* 4 + 3045t 3 + 3852t 2 + 2604f + 720 
U ~ 8z^ 2 (7i 2 + 24i+18) ' { ' 

At 7 + 60t 6 + 389t 5 + 1410f 4 + 3045t 3 + 3852t 2 + 2604^ + 720 
h ~ 8^(7i 2 + 24i+18) ■ ( ' 

The ratio of return times for N = 7 for three cases: 
for u — const. 

*l-W*i-h = (5184000t 16 + 90633600i 15 + 734038560* 14 + 3656624904i 13 + 12543798852i 12 + 31435490078< n + 
59579986661i 10 + 87223274254i 9 +99846813214t 8 + 89833419890i 7 +63379753809t 6 + 34652851894t 5 + 14319281196t 4 + 
4279417752t 3 + 859191840i 2 + 101520000i + 5184000)/(5184000i 10 + 60393600* 9 + 306948960i 8 + 896350104£ 7 + 
1664901648i 6 + 2053477662i 5 + 1700206878i 4 + 930252240i 3 + 320428800i 2 + 62380800t + 5184000), 

for Ui/v = 9/i 

^ L /t w = (576t 16 + 16800t 15 +229696t 14 +1956752i 13 +11645844t 12 +51472058t n +175326610i 10 +471411274i 9 + 
1015867913t 8 + 1768373403i 7 + 2486683328i 6 + 2797983376i 5 + 2465006400^ + 1636404624t 3 + 767257920t 2 + 
225504000i + 31104000)/(76032t 10 +1450368t 9 + 12336072t 8 + 61572600t 7 +199652130t 6 + 439389384t 5 + 664690536t 4 + 
682575840* 3 + 455457600< 2 + 178329600* + 31104000), 

and for a = a/i 2 

tL-n,/h-n = (27000t 16 + 1648350i 15 + 46021545i 14 + 781598610i 13 + 9060806565i 12 + 76286696592i u + 
484749056302t 10 + 2385421175676t 9 + 9253317988496i 8 + 28615082281632i 7 + 70836261328608i 6 + 
1396362454773m 5 + 215233793554176i 4 + 250177275371520i 3 + 205696375603200i 2 + 106205478912000t + 
25798901760000)/(1032264000t 10 + 28816738200i 9 + 362493838440i 8 + 2695926522960t 7 + 13095290178720t 6 + 
43360671643200t 5 + 99111840724224* 4 + 154547465674752t 3 + 157566016143360t 2 + 95025954816000t + 
25798901760000). 

The exact value of af for N = 10: 

E 2«2Q + ^43 + 2^ 4 4 + 2t; 45 + vajq + ^72 + 2t; 73 + v 7i + v 77 + 2v 79 + v g5 + v 96 

4 2(^20 + W43 + V44 + V 45 + V46 + v 72 + v 73 + v 74 + v 77 + V 7S + v 7g + v 95 + v 96 + v 98 + v 99 + 2v w6 ) ' 

A prime factorization of the common factor of probabilities of states of FRDA for N = 10 (presented in Table IXTV]) 
for rebound parameters u — 1 and v = 1: 

2 2 * 3 2 * 607 * 66617 * 35622218878023086289926346229321892786614786277415883867. 

The biggest prime has 56 digits. 
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TABLE XIV: All states for the size of the lattice JV = 10. The probabilities are for the case /i = 1 and v = 1. States which differ 
by mirror symmetry are indicated in the column with label 'sym'. Values in parenthesis indicate smaller than 10 translational 
symmetry. 
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